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\0 \ Abstract 



We present a formalism that allows the computation of the baryon asymmetry of the universe 
from first principles of statistical physics and quantum field theory that is applicable to certain types 



^ . of beyond the Standard Model physics (such as the neutrino Minimal Standard Model - z/MSM) 

and does not require the solution of Boltzmann or Kadanoff-Baym equations. The formalism 
^> ' works if a thermal bath of Standard Model particles is very weakly coupled to a new sector (sterile 

^ : 

CN . neutrinos in the z^MSM case) that is out-of-equilibrium. The key point that allows a computation 



without kinetic equations is that the number of sterile neutrinos produced during the relevant 
cosmological period remains small. In such a case, it is possible to expand the formal solution of 
the von Neumann equation perturbatively and obtain a master formula for the lepton asymmetry 
expressed in terms of non-equilibrium Wightman functions. The master formula neatly separates 



H ' CP-violating contributions from finite temperature correlation functions and satisfies all three 

Sakharov conditions. These correlation functions can then be evaluated perturbatively; the validity 
of the perturbative expansion depends on the parameters of the model considered. Here we choose 
a toy model (containing only two active and two sterile neutrinos) to illustrate the use of the 
formalism, but it could be applied to other models. 
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I. INTRODUCTION 



The baryon asymmetry of the universe (BAU) is a very important quantity in cosmol- 
ogy. It is quantified by the (dimensionless) baryon-to-photon ratio rj. This ratio has been 
measured very precisely by the WMAP collaboration [Lj and is given hj rj = {rii, — ni)/n^ = 
(6.1 ± 0.2) X 10~^°. Here (n^) is the number density of baryons (anti-baryons) and is 
the number density of photons. If this excess of baryons over anti-baryons is not simply an 
initial condition of our universe, then it is the goal of any particle physics models to explain 
this asymmetry. 

In principle, the Standard Model possesses all the necessary ingredients to produce an 
asymmetry (i.e. it fulfills the three Sakharov conditions baryon number violating 

processes are mediated by sphalerons, CP-violation is hidden in the Cabibbo Kobayashi 
Maskawa (CKM) matrix and a first order phase transition would provide the necessary out- 
of-equilibrium condition. On the other hand, it has been shown that in the Standard Model 
there is no first order electroweak phase transition for Higgs masses above 80 GeV (it is a 
smooth crossover) Since the lower bound on the Higgs mass from LEP is 114 GeV, the 
out-of-equilibrium condition necessary for baryogenesis is not satisfied and no asymmetry 
is produced. In extensions of the Standard Model, this conclusion may change; thus new 
physics beyond the Standard Model is necessary to explain the BAU. 

There exists many mechanisms/models that could explain the BAU. The most well-known 
are Grand Unified Theories (GUT) baryogenesis, electroweak baryogenesis in extensions of 
the Standard Model, leptogenesis and the Affleck-Dine mechanism. See Refs. fljioi for 
reviews. A common feature of these models is that they require heavy degrees of freedom 
that are hard to detect with present accelerators (the LHC may discover some of these 
particles in the near future). 

An alternative model is the neutrino Minimal Standard Model (i/MSM) ll|-|l9|: it is a 
minimal extension of the Standard Model with three sterile right-handed neutrinos with 
masses below the electroweak scale. This model has interesting features. For instance, it 
could explain simultaneously three shortcomings of the Standard Model (namely neutrino 
oscillations, dark matter and BAU). Also, since it contains no very heavy degree of freedom, 



it could in principle be tested experimentally with present facilities 
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2l|. 



For baryogenesis to occur some degrees of freedom must be out of thermal equilibrium 
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in order to fulfill the third Sakharov condition. The usual approach to baryon excess com- 
putations is to use Boltzmann equations. Various assumptions are used in the derivation 
of Boltzmann equations, one of them being that the coherence length of the processes in- 
volved must be much smaller than the mean free path of the particles. As long as these 
assumptions are satisfied, Boltzmann equations can describe systems that are arbitrarily 



out-of-equilibrium. But it is shown in Refs. 22|, |23| that coherence effects are important in 



baryogenesis, and thus a more refined quantum mechanical treatment is needed. 

Out-of-equilibrium quantum field theory is a notoriously difficult subject. Significant 



progress has been made recently in formal aspects (e.g. 2J-|26|) but applications to realistic 



baryogenesis computations are still lacking (although see the recent progress in Refs. 
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32|). 

Our utlimate goal is to compute the BAU in the phenomenologically interesting z/MSM. 
Previous calculations of the BAU in the i/MSM [12| show that the produced asymmetry 
is in the right range. The production of a baryon asymmetry in the z/MSM happens via 
coherent active-sterile neutrino oscillations and requires appropriate kinetic equations for its 
treatment. These calculations are performed using kinetic equations of the form 
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3^: 



= [H,p]-\{T,p]+'-{T,^-p}, (1) 

where p is the complete neutrino density matrix, H is the Hamiltonian and Tp and F are 
production and destruction rates, respectively. The approach of Ref. [12| based on Eq. ([1]) 
has several weak points. First, Eq. ([1]) relies on the usual assumptions of kinetic theory (with 
the additional assumption that the duration of a collision is small compared to the various 
oscillation times of p). Second, the calculations in Ref. 12| are done in the relaxation time 
approximation and it is assumed that the integrals in the collision term are dominated by 
0(T) momenta. This assumption might not be warranted; since this is basically an oscilla- 
tion problem with many particles, there are many timescales involved and all components 
might not relax in the same way. Another weak point of Eq. ([1]) is that it is not systemati- 
cally derived from first principles and thus there is no real control over the error. This has 
important consequences for phenomenology since "factors of a few" in the determination 
of the allowed parameter range of the z/MSM are crucial for the experimental searches of 
its particles. Thus in order to better constrain the model and study its phenomenological 
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implications, we need a first principles calculation based on quantum field theory. The first 
principles formalism presented here and the treatment based on Eq. ([1]) could also be com- 
pared in a region of parameter space where the two approaches should match in order to 
estimate the accuracy of kinetic theory approaches. 

In this paper, we present the first steps toward a computation of the baryon asymmetry 
from first principles of statistical physics and quantum field theory. We focus on a particular 
class of models similar to the z/MSM, namely a minimal extension of the Standard Model 
with an arbitrary number of sterile neutrinos. Our approach is sufficiently general that 
it could be applied to other models. The baryogenesis scenario studied here is similar in 
spirit to leptogenesis (with crucial differences). The key idea is that in some region of 
parameter space where the Yukawa couplings of the sterile neutrinos are small, it is possible 
to use conventional perturbation theory without having recourse to more sophisticated non- 
equilibrium tools (such as Kadanoff-Baym equations). For a similar treatment used in the 
context of sterile neutrino production, see Ref. [3]. 

The rest of the paper is organized as follows. In Sect. [Tll we present the model Lagrangian 
used in this study and outline our baryogenesis scenario. The core of the paper is Sect. IIIIl 
where we derive a perturbative formula that expresses the (lepton) asymmetry in terms 
of Wightman functions and discuss its range of validity. In Sect. llVt we test the formula 
derived in Sect. IIIIl by applying it to a toy model that can be solved both perturbatively 
and exactly. This computation also serves as an illustration of the inner workings of the 
perturbative formula. Some elements of this paper have already been published (in a different 
form) in the proceedings of Strong and Electroweak Matter 2008 35| . 

II. THEORETICAL BACKGROUND 
A. Model Lagrangian 

A generic Lagrangian containing A active neutrinos and B sterile (right handed) neutrinos 
can be written in the chiral basis as {a runs from 1 to A and / from 1 to B): 

Cab = L^tl^L^ + Nji^Ni - ^N^Nj - F^jL^Nj^ + h.c, (2) 

where L is the active lepton doublet, N the sterile neutrino singlet, A^"^ the charge conjugated 
A, $ the Higgs doublet, p the Higgs field, M/j are Majorana masses and F^i are Yukawa 



TABLE I: 



A B Mi J Yuk. + Mix. Phases Total 

3 3 3 9 6 18 

3 2 2 6 3 11 

3 11 3 4 

2 2 2 4 2 8 

2 11 2 3 

1 2 2 2 1 5 



couplings. Here sterile means singlet under the Standard Model gauge group. Note that for 
the purpose of deriving a master formula for the asymmetry, only active- sterile transitions 
are necessary and other Standard Model particles are left out here. The case A = 3 and 
5 = 3 (and including Standard Model particles) corresponds to the z/MSM. Note that 
in the phenomenologically relevant region of the z/MSM parameter space where the right 
abundance of dark matter is obtained, the "dark matter sterile neutrino" has a tiny Yukawa 
coup.., e.s.„t.ny decouple, Q. TKus the .MSM has effectively tK..ee active and 
two sterile neutrinos. 

We can check that the Lagrangian ([2]) satisfies all three Sakarov's conditions for baryogen- 
esis. First, the presence of the Majorana mass term breaks lepton number conservation (and 
thus baryon number conservation via sphaleron processes). Second, if Faj is complex, then 
there can be CP violation in the system (similar to the case of the CKM matrix). Table (jl]) 
shows the counting of physical parameters (masses, Yukawa couplings, mixing angles and 
comlex phase) in (|2]) for various numbers of active and sterile neutrinos. Note that complex 
phases can only be present for B > 2. CP violation is thus possible in the z/MSM. Third, 
the out-of-equilibrium condition is provided here by the expansion of the universe. If the 
rate of sterile neutrino production Fsteriie is less than the rate of expansion of the universe 
H, then sterile neutrinos do not interact enough to preserve equilibrium. This ratio can be 
estimated as follows: 

-Tsterile / 



where we used T = Tgph ~ 100 GeV (temperature at which sphalerons become inefficient) 
and / = 10~^ (a possible Yukawa coupling in the z/MSM). Since the ratio is inversely propor- 



tional to the temperature, the sterile neutrinos are out-of-equihbrium for all temperatures 
of interest. Note that this ratio depends on / and that the Yukawa couplings in the z/MSM 
are constrained by observations 17|. The value used in the above estimate is for a typi- 
cal choice of parameters; there is a choice of coupling when sterile neutrinos equilibrate at 
T > 100 GeV. In the latter case our formalism does not work and kinetic equations (or more 
sophisticated non-equihbrium quantum field theory tools) must be used. 



B. Outline of the Baryogenesis Scenario 

The scenario for baryon asymmetry generation in sterile neutrino oscillations was pro- 



posed in Ref. |3J] and developed in Refs 



12 



17|. Just after reheating (tj = 0), we have 



a thermal distribution of Standard Model particles and we assume that there is no sterile 
neutrino initially. This is true if there is no source of sterile neutrino during inflation, such 
as in recent models where the inflaton is played by the Higgs 15|, |36|. There are no ex- 
perimental data that could constrain the initial conditions at the moment. The final result 
for the asymmetry of course depends on the initial conditions. Note that the formalism 
presented in Sect. IIIII is valid for any initial distribution of sterile neutrinos. 

During the following cosmological evolution, the three Sakharov conditions are satisfied 
and a lepton asymmetry is produced via coherent and resonant oscillations of active-sterile 
neutrinos. The necessary resonance condition for a sufficient asymmetry generation is that 
two of the sterile neutrinos should be nearly degenerate in mass. This lepton asymmetry 
is converted into a baryon asymmetry via sphaleron processes. Since sphalerons become 
inefficient at a temperature around Tgph ~ 100 GeV, the conversion process stops at around 
t/ = tsph = 10i' GeV-\ 

This scenario is similar in spirit to "usual" or thermal leptogenesis, although the physics 
is different. In the case of thermal leptogenesis, the parameters of the model are different 
(large Yukawa couplings and large masses for the sterile neutrinos) and the sterile neutrinos 
are initially in thermal equilibrium. It is after this initial period of thermal equilibrium 
that they go out-of-equilibrium and decay, producing a lepton asymmetry. This period of 
out-of-equilibrium is relatively short. In the case of the z/MSM, sterile neutrinos are not in 
equilibrium initially and typically stay out-of-equilibrium (due to their very small couplings) 
until sphaleron processes become inefficient. The exact time at which sterile neutrinos reach 



6 



thermal equilibrium depends on the parameters of the model (see Sect. IIII Bp . 



III. PERTURBATIVE FORMULA FOR THE (LEPTON) ASYMMETRY 
A. Derivation of the Master Formula 

The total lepton asymmetry per unit volume is given by the difference between the 
average number density of leptons minus the average number density of anti-leptons. The 
asymmetry at time t for a particular flavor a of active neutrino z/^ is: 

A„(t,f) = Tr [p(t)z/t(t,x)z/„(t,x)] . (4) 

All fields in Eq. (jl]) are in the interaction picture. Here p(t) is the appropriate density 
operator. The initial condition for the density matrix is expressed as: 

p(0) = ps®PsM, (5) 

where ps is the out-of-equilibrium density matrix for the sterile neutrinos and psM = e~^^^/"^ 
is the usual equilibrium density operator for Standard Model particles. The precise form of 
p'g is irrelevant for the derivation of the master formula, but the final value of the asymmetry 
is dependent on it. For the case of the baryogenesis scenario outlined in Sect. Ill B] the density 
matrix for sterile neutrinos is ps = |0)(0| at t = 0, where |0) is a sterile neutrino vacuum 
state. 

Equations (jlj) and ([5]) constitutes the system that we want to solve (for any particular 
flavor). In equilibrium, the density operator is an exponential and is time- independent; the 
machinery to solve such problems is very well developed. In the present case, Eq. ([3]) shows 
that the system is out-of-equilibrium, meaning that the density operator is time- dependent 
and that usual equilibrium techniques fail. A simple way to see why such problems are hard 
is that propagators at finite temperature depend on distribution functions. The building 
blocks of perturbation theory are propagators and vertices. In out-of-equilibrium situations, 
distribution functions are time- dependent, implying that propagators change with time in 
a non-trivial way. Thus usual perturbation theory fails and resummation techniques must 
be used [26|, |37|, [SSj , unless changes in the propagators are small on timescales relevant for 
the problem at hand. In the following, we treat Eqs. (Il])-([5]) perturbatively and show in the 
next section under what condition perturbation theory is valid. 



The first step is to find tlie time evolution of tlie density operator. It is given by tlie von 
Neumann equation: 

.dpit) 



dt 



Hi{t),m 



(6) 



wliere H{t) = Ho+Hj{t) is the interaction Hamiltonian corresponding to the Lagrangian (E]). 
All operators in Eq. are in the interaction picture. Note that we work in fiat spacetime 
here; effects due to an expanding background can be incorporated using the usual procedure 
of expressing the equations in conformal time and going back to physical time at the end of 
the computation. The rest of the derivation does not depend on the background. 
We can find an iterative solution to Eq. (j6]): 



m = p(0)-z / dt' Hiit'),p{0) 



- t dt' t dt" 

Jo Jo 



Hjit'), Hjit"),p{t") +0{H!).{7) 



Note that the expansion parameter in Eq. ([7]) is ^^Hjt" . Thus the iterative solution ([71) is 
also a perturbative solution if the criterion ''Hjt <^ 1" is satisfied. We use quotation marks 
here to indicate that this is not a precise criterion since Hj is an operator and t is a number; 
we present a more precise criterion in the next section. 

Substituting the iterative solution (I7l) back into Eq. (jll), we schematically get: 



A,(t',x) = TT[{m + p'^'\t')+p('\t') + ...)ui{t',x)iy^{t',x) 



(8) 



The first contribution is just an irrelevant infinite constant independent of the temperature; 
it is an artefact of using the lepton current to compute the asymmetry and we discard 
it in the following. The second contribution (and all contributions p*^"^ with n odd) are 
automatically zero because they contain an odd number of creation/annihilation operators. 
The first non-trivial contribution is thus (we only keep terms up to 0{Hj) in the following): 

rt' 



A^{t',x) 



dt / dti Tr 
Jo 



iHj{t),[Hj{t^),p{0)]])uiit',xM',x 



(9) 



To make progress, the interaction Hamiltonian need to be specified. A convenient form for 
our purposes is to decompose it as: 

Hj{t) = J d'x (vp{t, x)Jp{t, x) + J^{t, x)v^{t, x)) , (10) 

where J contains all the information about the interaction except the field associated to 
the desired asymmetry. The exact form of J is not important for the rest of the derivation. 
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Note that J contains the out-of-equihbrium sterile neutrino field; in the case of the z/MSM, 
it also contains the Higgs field and the Yukawa couplings. 

Inserting this form for the interaction Hamiltonian into Eq. ^ and using Wick's theorem, 
we obtain: 

-{J^Jp){y-yyi){yayp) + {J'^Jp){yiyp){h'-yi'a) +0{Hf), (11) 

where we have used a condensed notation J/j = Ji^ait^y) with the first and second indices 
being fiavor and spinor indices, respectively; a always refers to the external fiavor index and 
spacetime coordinates (t', x), while /3, 7, ... correspond to internal fiavor/spinor indices and to 
spacetime coordinates (t,y), (ti,yi), .... The equilibrium correlators {uu) are the usual active 
neutrino propagators. The correlators {J J) contain the out-of-equilibrium sterile neutrino 
fields and potentially other fields $ that are in equilibrium; they can be decomposed as 
{J J) = Tr \ps PsmJJ] ~ Tr \psNN] O Tr [psM$5> 



Note that for the scenario outlined 



in Sect. Ill Bl we have ps = |0)(0| and thus zero temperature propagators can be used to 
describe the evolution of the sterile neutrinos. This replacement is allowed if perturbation 
theory is valid; we come back to this point in the next section. 

It is possible to simplify Eq. ( ITTl) further by making an assumption. First note that for 
any propagator we have {il){x)il){y)) ~ {ijj{x)ilj{y)), where the ~ means everything is equal 
except for their Dirac structures which are complex conjugate of each other. Looking at 
Eq. ffTTl) . we see that the first and second terms and the third and fourth terms are the same 
except that their coupling constants and their Dirac structures are complex conjugate of 
each other. In the following, we assume that the Dirac structure computation for each term 
gives a real scalar; thus for all practical purposes, the Dirac structure is irrelevant when 
comparing these pairs of terms. This assumption may be true in general, but its validity 
must be checked explicitly for each interaction Hamiltonian. 

Taking into account the previous assumption and defining Jg = Fjsjjj, the asymmetry 
at second order in the interaction Hamiltonian finally becomes: 

A„(t',f) = -J^ dtj^dh J d'y J d''y,2Re[2tlm{FpiF;j){JiJj){ulu.,){u^up) 

-2iIm{F;jF^j){JjJj){iy,iyi){u^iy^) 



dtj^dh I d^y J dSi 4 \lm{F^iF;j) \m{{Ji'j j){viv^){v^vp)) 

-\m{F;jFp,) \m{OjJ,){v,vi){v^vp))\ . (12) 

This last equation, or more specifically its 0{Hf) version (see the discussion below), are 
the main results of this paper. It expresses the lepton asymmetry for a particular flavor in 
terms of Wigthman correlators. If perturbation theory is valid, then these correlators can 
be computed using conventional tools. If the result is zero, then higher order contributions 
in Hj must be computed. 

Formula f lT2|) is quite general, applicable in principle to leptogenesis-type models. The 
validity of perturbation theory is the only assumption that enters into its derivation (in 
addition to the minor assumption about the Dirac structure). 

The last formula satisfies the third Sakharov condition. As mentioned previously, the J 
operator contains the sterile neutrino operator, and the sterile neutrinos are not in thermal 
equilibrium. Assuming that the sterile neutrinos are in equilibrium, then the J operator 



would obey the Kubo-Martin-Schwinger condition [39|,|40| {Ja{ti)Jb{t2)) = {Jb{t2)Ja{ti+'i/^)) 
[P = 1/T is the inverse temperature). Using the Kubo-Martin-Schwinger condition and 
time translation invariance, it can be shown that the asymmetry is automatically zero in 
equilibrium. 

Formula ( !T2|) also neatly separates CP violating effects (contained in the imaginary part 
of the couplings) and dynamical effects. It is thus easy to see that it satisfies the second 
Sakharov condition: if there is no CP violating phase, then the Yukawa couplings F's are 
real and the formula gives automatically zero. Furthermore, using the fact that propagators 
are diagonal in flavor space ((z/aZ//?) oc 6a[s and {JiJj) oc (5/j), we get that Im [PpiFpj^ = 
and that the formula for the asymmetry is automatically zero at 0{H]). This is another 
way of saying that a lepton asymmetry is a quantum effect (generally coming from the 
interference of a tree level diagram and a loop diagram). 

The first non-trivial order in the expansion of the density operator is thus 0{Hf). Ex- 
panding the iterative solution ([7]) up to 0{Hf) and repeating the same procedure, we finally 
obtain: 

A^{t',x) = - f dt f dti r dt2 r dh f d^y I Ai / ^2 / As 

JQ Jo Jo J J J J 

x4 [lm(F^,F^jF;^F;^) Im[( J, JjJ^^Ji) 
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+\m{FpiF;jFsKF:L) lm[{JjJjJKJL) 
+lm{FpiF;jF*,j,F,L) lm[{JjJjJKJL) 
+Im(F,lF^,F^jF,V) lm[{JLJiJjJK) 
+lmiF:^FpiF;jFsK) \m[{J Jfj jJk) 

{-(^e^Di^-yi^s) - {^i^j) {^ei^s)) {{^p^a) + {^aVp)) 

+lm{F,LFpiF;jF*,K) \m[{jJiJjJK) 
+lm{Fi^FpjF^jF:^) ImlCjKJiJjJL) 
+lm{FiKF[sjF;jF,L) \m.[{J kJiJ jJl) 
+\m{FsKFpiF;jF:^) lm[{JKJiJjJL) 
+lm{F:^Fi^FisjF^j) lm^[{J J kJiJj) 
+lm{F,LFiKFpjF;j) ImlijjKJiJj) 
+lm{F:^FsKFisjF;j) lm^[{J JkJiJ j) 



(13) 



Formula f|T3|) allows the computation of the BAU up to 0{f^) and is exact in all other 
couplings of the model. Since perturbation theory is assumed to be valid, it is possible to 
resolve the 4-point functions (JiJjJrJl) into products of 2-point functions using Wick's 
theorem. Thus the only inputs in Eq. (IT^ are the -F^/'s and the various propagators of the 
model. The resulting terms can be interpreted as Feynman diagrams. 
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N V N V V \ N V , N V vN v N 

H 



FIG. 1: Typical Feynman diagrams appearing in the perturbative formula (jl3p . The solid lines 
correspond to active v and sterile neutrino propagators and the dashed lines correspond to Higgs 
H propagators. 

The creation of the asymmetry takes place in the early universe where a thermal bath 
of Standard Model particles is present. It is well known that resummed propagators (with 
thermal masses and dampings) must be used in order to obtain correct results. For instance, 
interactions of particles with the medium may open new decay channels that are otherwise 



kinematically forbidden (see for example 4l|]). The inclusion of damping is also important to 



deal with so-called secular terms or "pinch" singularities (see Sect. IIIIBj) . Hard thermal loop 



42| 



resummations are moreover necessary to obtain gauge invariant results in some cases 

To end this section, we discuss qualitatively the case of the z/MSM in order to illustrate 
the use of the perturbative formula ( |T3l) (for a similar discussion in a kinetic theory setup, 
see [l7|). We restrict ourselves to the symmetric phase here. The typical Feynman 
diagrams appearing in the perturbative formula are given in Fig. ([1]). In order to compute 
the asymmetry, we need the (Wightman) active/sterile neutrino and Higgs propagators. 
The exact form of these propagators depend on the self-energies that are resummed into 
them. For sterile neutrinos, no resummation is necessary since the propagators are at zero 
temperature and the self-energy corrections are small (i.e. 0{p)). For active neutrinos, it is 
necessary to include W, Z boson and charged lepton loops. The absence of charged leptons 
would imply a new symmetry for the z/MSM, allowing the removal of CP phases in Fai by 
rotating active neutrino fields and leading to a vanishing asymmetry [51]. For Higgs bosons, 
the dominant contribution to the self-energy comes from top quark loops of size 0(m^T^/f ^); 
this is parametrically large compared to W,Z boson loops of size 0{m\^ ^h"^) Higgs 
loops of size 0{m\j/v'^) (all the estimates are for large temperatures). The self-energies that 
need to be resummed in the propagators are illustrated in Fig. ([2]). 

As can be seen from the above qualitative discussion, the computation of the BAU in the 
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FIG. 2: Self-energies that are resummed in the active neutrino v and Higgs H propagators. 

z/MSM is not a simple task. For this reason and as a warm up, we analyze a simpler toy 
model involving only Yukawa interactions with sterile neutrinos and no other interaction. 
The full computation of the BAU in the z/MSM will be the subject of a separate publication. 



B. Validity of Perturbation Theory 



For the iterative solution ([7]) to be also a perturbative solution, the criterion "iZ/t <^ 
1" needs to be satisfied. Said differently, it means that perturbation theory is bound to 

break down for sufficiently long times; these are the infamous secular terms that plague 

I ' 

non-equilibrium quantum field theory (e.g. [26]). These secular terms also appear in a 
different form in quantum field theory computations using the real-time formalism of finite 
temperature field theory 43|. In the real-time formalism (where time is taken to go from 
minus infinity to plus infinity and back), ill-defined products of delta functions (or "pinch" 
singularities) arise naturally [44, 



45] and blow up with the time- volume 46|]. Fortunately 



;hese pinch singularities cancel in equilibrium due to the Kubo-Martin-Schwinger condition 



471]. Since the Kubo-Martin-Schwinger condition is not valid out-of-equilibrium, the pinch 



singularity; 



be used 



prob 



em remains for non-thermal systems and more sophisticated methods must 



38]. 



In the present paper, we use the von Neumann equation and time is finite: if the final 
time is taken to be small, then secular terms should also stay small (see j^sl for a similar 
argument phrased in terms of pinch singularities). We can estimate the size of the secular 
terms in the following way. Roughly speaking, "ff/t" gives the number of sterile neutrinos 
produced per unit time {Hi) times the time (t). The total number of sterile neutrinos 
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produced is given by: 

y-, / ,/N ,,/ power laws i sterile 



'^sterile ~ / rsterile(^ )dt 

Jo 



(14) 

^sph 



H 

where in the last step we use the fact that both the rates of sterile neutrino production and of 
the universe's expansion are power laws. Thus the size of secular terms (or equivalently the 
total number of sterile neutrinos produced) is given by the out-of-equilibrium criterion ([3]) 
evaluated at the time when sphalerons become inefficient. For T = Tgph = 100 GeV and 
/ = 10~^, we get nsteriie ~ 0.1 ^ 1 and thus perturbation theory is justified in this part 
of parameter space. It also shows that the distribution function for the sterile neutrinos 
does not evolve much during the relevant period for baryogenesis and justifies our use of 
zero temperature propagators to describe their evolution (in the case of the baryogenesis 
scenario outlined in Sect. Ill Bl) . 



IV. APPLICATION TO A TOY MODEL 



Before using the formula for the asymmetry ( fT3ll on a realistic but more complicated 
model (such as the i/MSM), we would like to test it on a simpler (but unrealistic) model and 
verify the range of validity of perturbation theory. We thus make the following simplifications 
to the model Lagrangian ([2]). First, we neglect interactions (damping) with Standard Model 
particles. This implies that the total active lepton asymmetry is zero (see the discussion 
below Eq. f[T^ ) and only individual fiavor asymmetries are non-zero. Second, we assume that 
the Higgs is non-dynamical, making the Lagrangian quadratic in the fields. This simplified 
model thus describes a bunch of non-interacting harmonic oscillators and no thermalization 
is possible. Third, we reduce the number of free parameters by reducing the number of 
active/sterile neutrinos. Looking at Table ([T]), the minimal model containing only one CP 
violating phase has one active and two sterile neutrinos. For reasons that will become clear 
in the next section, the case with one active neutrino and two sterile neutrinos trivially gives 
a vanishing asymmetry. We thus opt for the next-to-minimal model (i.e. two active and two 
sterile neutrinos) to test our formula. The Lagrangian for such a toy model is: 

-Mi)l^ia'^i)2R ^(^?R^(^^^iR + ^L^(^^V'2i?) 

-fvi>lL^iR - efve''^i>lL'ilj2R - 6fviljlLi>2R + h.c, (15) 
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where we take the Higgs field to be a constant v, ipi,2R are the sterile neutrino right-handed 
Weyl fields, V'3,4-L are the active neutrino left-handed Weyl fields, M is the common mass of 
the sterile neutrinos, AM is the sterile neutrino mass difference, / and e and 6 are Yukawa 
couplings and 1] is a CP violating phase. 

This toy model has five real parameters and one CP violating phase [52] (compare this 
to the 8 real parameters and 3 phases of the i/MSM including constraints from dark matter 
abundance). The physics of leptogenesis in this toy model is thus simpler. For instance, if 
?7 or e or AM are sent to zero, then the Lagrangian does not contain any CP or baryon 
number violating terms anymore and the asymmetry vanishes (if 6 is zero then we come 
back to the one active neutrino case and the asymmetry is also zero). These features should 
also appear in the solution. 

We can solve this model in two ways. Since the Higgs takes its expectation value, the 
Lagrangian is quadratic in the fields and it is possible to solve the system exactly. If the 
Yukawa coupling / is small, then the system can also be solved perturbatively using the 
master formula for the lepton asymmetry presented in Sect. IIII A[ The exact and pertur- 
bative results can be compared and should match in some time interval. This is what we 
present in the next sections. 



A. Exact Solution 

As in the perturbative case, the lepton asymmetry (for the active fiavor a = 3, 4) is given 
by the average lepton current density: 

A,(t,f) = Tr [p(0)7/'ii(t,f)z^„L(t,f)] , (16) 

where all fields are in the Heisenberg picture. The calculation of the asymmetry is now purely 



quantum mechanical and very similar to neutrino oscillation computations (e.g. j49|). For 
times t < 0, there is no interaction and there is no sterile neutrino. At t = 0, the interaction 
is adiabatically switched on and sterile neutrino fields are initially in fiavor eigenstates (in 
which the thermal averages are defined). This last point is very important because otherwise 
the problem is trivial. Indeed, since the Lagrangian f|T5l) can be diagonalized such that all 
four masses are real (see below), then it means that all CP violating phases can be absorbed 
and there is no asymmetry if the thermal averages are defined in the mass basis. In the 
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case where the thermal averages are defined in the initial fiavor basis, the asymmetry is in 
principle non-zero. 

The goal is to compute the asymmetry at some time t > 0, hence the fields ^Jsiityx) 
must be time evolved from to some time t. Time evolution in quantum mechanics is most 
conveniently done using energy eigenstates. In order to obtain the energy eigenstates, we 
need the exact masses. The mass matrix Ai in ffTSl) is: 



/ 



M 



\ 



V 



(17) 



AM M fv 
M AM efve-''^ 6fv 
fv efve-'^i 
Sfv 

The above symmetric complex mass matrix can be diagonalized using Takagi's factorization 
Ai = UVU^ where U is unitary (e.g. jsO]). The results for the masses are: 



I 



fhi 


= M + 


rfi2 


= M + 




_ Pv' 


fhs 


M 




= i 



2M 
/V 
2M 



:i + e' + 5' + 2p'), 
'l + e' + 6'-2p'), 



^' , ) . (18) 

where we used the parametrization AM = ( '^ ^m"^*^"* ) ^ defined pe*^ = ^ee*'' + (e + S)k. 
Flavor and energy eigenstates are related as ipR = Ui/jr and ipL = i^*i^L where ipR = 
ii^iR iJ2R icF^i^li icr'^i^lLf and V'l = (-^o'^^ir - «'7^V^2fl ^3l)^- In the basis where 
the mass matrix is diagonal, the neutrino fields are given by: 
d?p 1 



(2vr) 



2EI h±i 



(19) 

where h is the helicity and the X(/i)(p)'s are two-component helicity eigenstate spinors. The 
exact energies are given by = ^Jlp]"^ + and the exact masses rfij's are given in Eq. f|T8|) . 

After t = 0, the initial flavor eigenstates become linear combinations of energy eigenstates 
and start to oscillate, with each energy eigenstate oscillating with a different frequency 
depending on the exact masses. At time t, we have: 



A„(t,f) = Tr p{0)[i>i{t,x)U^UU*Mt,x)]a 



(20) 
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In the last expression, the neutrino fields are expressed in terms of creation/annihilation 
operators dp^h,i that create/annihilate quanta corresponding to states with a certain mass. 
To do the thermal averages, it is necessary to re-express these "mass" creation/annihilation 
operators as linear combinations of "fiavor" creation/annihilation operators ap h,i- This is 
done using Bogoliubov-type transformations. The appropriate transformations are: 



E,+ \p\h 



p,n,t 



2E, 



j=l,2 



h,j\/Ej + \p\h- h(f){p, h)alp j^ j^Ej - \p\h 



fc=3,4 



^JE^-\v\h 



2E, 



i=i,2 



2E, 



JEj - \p\h + h(j){p, h)al ,^JEj + \p\h 



A;=3,4 



(21) 



where V is the matrix that diagonalizes the free part (i.e. / = 0) of the Lagrangian ([T5j) . 
Ei = \J\p\'^ + mf are the eigenenergies corresponding to the masses of the free Lagrangian 
and (j){'p,h) is some phase that satisfies (pi^—p.—h) = (f)*{p,h) and (f){p,—h) = —(f)*{p,h) 
491]. One can verify that these transformations preserve the canonical anti-commutation 
relations. 

Inserting the neutrino field operator f|T9|) into Eq. fl20|) and using the relations between 
the two sets of creation/annihilation operators (pTj) . the asymmetry becomes (after some 
algebra) : 



A„(t,f) 



2+A 2+A 

-EE 

i<j=l K=3 



(fp 
(2^ 



2b1 



^n{E,)\m{UM:,Um 



KJ 



[{E, + Ej) sin {E, - Ej)t - (Ei - Ej) sin (E^ + Ej)t)]{22) 



This is the final result for the exact lepton asymmetry. The above formula is true for two 
sterile neutrinos and any number of active neutrino fiavors A. The solution contains a CP 
structure part and a dynamical part consisting of two oscillating functions (one with a large 
amplitude and small frequency and one with a small amplitude and large frequency) for each 
value of We note that the CP structure part Im (UaMaj^Ki^Kj^ is similar to the one 



obtained in Ref. 121]. We also immediately see that the total active lepton number is zero, 
i.e. Ea Aa(t) oc Ea E-</=i Eltj Im (UaiU*jU*iU^j) = 0. In particular, we get A«(t) = for 
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only one flavor of active neutrinos; this explain our study of the two active and two sterile 
neutrino case. 



B. Perturbative Solution 

If the Yukawa coupling / is small, then the active-sterile neutrino interactions are 
small and can be treated as a perturbation. Diagonalizing the quadratic part of the La- 
grangian (fTS!) with the transformations ipiR {jiir + irj2R) / \/2 and iI)2R — ?■ (?7i_r, — «?72r) / 
(the active neutrino fields stay the same V'3,4L ~^ Vsal)^ obtain: 

Aoy = ViRi(^^df,r]iR + riljii(T^d^r]2R + vli^^^d^i^^iL + vlLi^''d^,r]4L 
(M + AM) ^ 2 (M-AM) r 2 
2 ViR 2 V2R 

fv f . . , e/fe*'' + , . , Sfv f , . X , /ooN 

-^VkiViR + iV2R) ^j^vLiViR - ^V2R) - ^VidViR - iV2R) + h.c, (23) 

The perturbative masses are mi^2 = M± AM for the sterile neutrinos and zero for the active 
neutrinos. The interaction Hamiltonian can be obtained from the above Lagrangian: 

Hi = vIlJs + 4V3L + vliJi + 4viL, (24) 
where the J operators are given by: 

Js = ^{l + ee'^)r],n+'-^{l-ee''')r]2R = Fs.ViR + F32V2R, 
Sfv i6fv 

Ja = ~^'^^R ~ = F41VIR + F42V2R- (25) 

Substituting the operators ( 125!) in Eq. (fT3|) . it is a straightforward but tedious exercice to 
compute the asymmetry. In the following we present the computation of (a part of) the 
first term in Eq. (fT3l) : the others are done in a similar way. Note that the assumption about 
the Dirac structure is verified here (see the discussion following Eq. flTT]) ). The first term of 
Eq. is: 

A„(t',x)i = - / dt f dti f ' dt2 f ' dt3 f d^y I d^yi f d^y2 f d^y^ 
Jo Jo Jo Jo J J J J 



-Aim (FaiF^jF*j^F*j^ ) im {{0\Jai{y)Myi)JlKiy2)JUy3m 

ivLe ix)VaLc (Z/2) ) {v\Lbiyi)V-,Ld{y3)) ( {vLa iy)VaLe (x) ) + {T]^Le {x)r]l^^ (?/) ) ) 

+4Im (FaiF^jF*j^F*j^ ) Im {{0\Jai{y)Jbj{yi)JlKiy2)JUy3m 

{vile {x)T]aLd (ys) ) {v\Lbiyi)VlLciy2)) ( {vLa {y)VaLe (x) ) + {r]aLe (a;)r/La (?/) ) ) 



(26) 
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where a, b, e are spinor indices, a, 7 refer to active neutrino flavors and /, J, K, N refer to 
sterile neutrino flavors. The case where a = 7 corresponds to one active neutrino flavor only 
and gives zero. For definiteness we thus consider a = 3 and 7 = 4 in the following. Keeping 
only non-zero contributions coming from the imaginary part of the coupling constants, we 
obtain: 

Aa(t',x)i = - / dt f dti r dt2 r dh I d^y f d^yi f d% f d^y^ 
Jo Jo Jo Jo J J J J 

-4im (F3ii^4ii^32^42) Im {{0\viRa{y)viRbiyi)vlRc{y2)vlRdiy3)\0) 

{vlLeix)V3Lc{y2)){vlLbiyi)V4Ld{y3)){{vlLaiy)V3Le{x)) + {mLe{x)7]lL^{y)))) 

-4im {f,,f,2f;,f:,) Im {{0\viRa{y)v2Rb{yi)vlRciy2)vlRdiy3m 

ivLe ix)r]3Lc{y2)) {vlLbiyi)ViLd{y3)) ( {vLa (l/)^3Le (x) ) + {rf^Le {x)t]1^^ (?/) ) ) 

+4lm{FsiF^,F;,F:,) Im {{0\viRa{y)ViRb{yi)vlRc{y2)vlRdiy3m 

(vile {x)r]sLd (Z/3) ) {vlLb{yi)V4Lc{y2)) ( {vLa (l/)^3Le (x) ) + {r]3Le (a;)r/La (l/) ) ) 

+4im {FsAf;,f:,) Im {{0\viRa{y)mRb{yi)vlRc{y2)vlRd{y3m 

{vlLeix)V3Ld{y3)){vlLbiyi)V4Lc{y4)){{vlLaiy)V3Le{x)) + {mLe{x)7]lL^{y))) 

+ (same terms with IR ^ 2R and F^i Fa2) ■ (27) 

Assuming the validity of perturbation theory, all the neutrino fields are free fields and we 
can use Wick's theorem to decompose the above 4-point functions into products of 2-point 
functions. Since the system is translationally invariant, we can do the Fourier transform 
over space. Concentrating on the first of the eight terms in Eq. fl27|) . we have: 

Aa(t',x)i,i = f dt f dh r dt2 r dh I -^^Alm{FsiF^iF;^F:,) 
Jo Jo Jo Jo J (zTrj'^ 

^^[{vlLeV3Lc){p) {0\vlRcVlRd\^)i.P) {vllbViLd) {-p) 

{^WRarilRbm-mrilLaV3Lem + (r/3Le^L) ) (P)] • (28) 

To make progress, we need the form of the propagators for massless left-handed active neu- 
trinos and for massive right-handed sterile neutrinos. These are obtained from the expansion 
of the fields in terms of creation/annihilation operators (cf. Eq. (fT9|) ). The propagators are: 

{vdhWdtMp) = (l-^(^p))X(-i)(p)x}-i)(p1e-^''(*-*^) 
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^-^p h 



-iEp{ti-t2) 



+n{E,)hMxi-h){-p)X{h){-p)e'''-^'^-'^^ 



-t2) 



-'P h 



+n{E,)hMx\^^{-p)x\_^^{-p)e 



(29) 



Inserting these propagators in Eq. f l28|) . we obtain: 



A(t',f)i,i 



dt / dti / c/t2 / (^^3 
Jo Jo JO 



(27r)3 4Im(F3iF4iF3*2F;2, 



Im 



mim2 
AE1E2 



((l-n(E3))(l-n(E4))< 



,i{-Ei-E3)t i{Ei-E4)tlJ{-E2+E3)t2j{E2+E4)t3 



+n{E3)n{Ei)e 



i{-Ei+E3)tJ{Ei+E4)tiJ{-E2-E3)t2j{E2-E4)t3 



(30) 



The integrals over time are trivially done; the result is: 

d^P 771x1712 



(27r)3 AE1E2 



t' 



l-n{Es)){l-7i{E,)) 

sin(Ei + ^3)^' 



{El + E3){E2 + E^){E:, + E^) {Ei + ^3)^(^1 - E^){Ei -E2 + E3- E^) 
sin(E3 + Ei)t' sin(E2 + Ei)t' 



+ 



(E3 + ^4)2(^1 - E^){E2 - Es) {E2 + E,Y{E2 - ^3)(^l -E2 + E3- E,) 



+n(£'3)n(i?4) (same terms with i^'3^4 — —-£^3,4) 



(31) 



The other terms in Eq. (1131) can be done in a similar way and they all have similar structures 
(with different arrangements of energies in the sines and denominators) . There are 12 x 4 x 
2 = 96 terms similar to Eq. fl3T|) in the final result for the perturbative lepton asymmetry. 
Because of its size (there are no obvious simplifications) and since it is not very instructive, 
we do not write the full expression of the perturbative lepton asymmetry here. 

Even at the level of indivitual terms we see that Eq. ( |3T1) has the features expected from 
the first two Sakharov conditions. First it is clear that Eq. ( 13T|) is zero when = 0, e = or 
6 = (i.e. no CP violation). It also vanishes when the sterile neutrino masses are degenerate 
(i.e. no lepton number violation). To see that, note that each terms in Eq. (I3T!) is paired up 
with a similar term with IR 4-> 2R and F^i o Fa2 (cf. Eq. fl27|) ). The asymmetry is thus 
proportional to Ei — E2 and vanishes when the two energies are equal (i.e. when AM = 0). 
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C. Comparison Between Exact and Perturbative Solutions 

The exact and perturbative computations of the asymmetry start with the same initial 
condition and the dynamics is dictated by the same Lagrangian. They should therefore give 
the same result up to some time where secular terms become important. To estimate this 
time, it is not possible to use the considerations of Sect. IIIIBI because there is no interaction 
(thus no thermalization) and spacetime is not expanding. We use instead the following 
argument. 

The toy model considered in Sect. IIVI is quadratic in the field, thus the asymmetry 
production should be oscillatory. Looking at the exact ( 12^ and perturbative (ISTi) solutions, 
we note that both solutions are sums of oscillatory functions with different frequencies. The 
"exact" and "perturbative" frequencies are given by ujij^±. = {Ei±Ej) and ujij^±. = {Ei±Ej). 
Since the masses are different in the exact and perturbative cases, the frequencies are also 
different. This implies that the two solutions develop a phase difference over time. This 
phase difference is secular. Thus even in a non-interacting theory, secular terms are present 
because of the building up of phase diffence between solutions. 

The exact and perturbative solutions should agree when this phase difference is small. 
We estimate this phase difference in the following. The frequencies can be approximated as 
^ij,± ~ i 1) + (^i i ^j)/2|p| (we assume that m < \p\ here). The criterion for the 
smallness of secular terms is thus: 



Thus the set of masses that produces the largest frequency difference gives the lowest time 
at which the two solutions should differ from each other. 

In the following we compare the exact result ( 122|) and the complete perturbative result (c.f. 
Eq. ( l3T|) ) numerically. We plot both results for representative time intervals. For simplicity 
we plot the asymmetry per unit phase space volume and take a typical momentum p ~ T. 
The values of the parameters used are / = 3 x 10^^, v = 174 GeV, e = 0.8, 6 = 0.3, r] = 0.6, 
K = 1.1, M = 1.5 GeV, T = 100 GeV and / = 10^^. The results are shown in Fig. (jS]). 

We can estimate the time tg at which secular terms become important using the crite- 
rion fl32|) and the expression for the masses f|T8|) . For the parameters chosen here, the largest 



Uij^± — Uij^± \t =^ t <^ 



m 



(32) 



(m; 



;f ± mj) - {mf ± m|) 
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FIG. 3: Plots of active lepton asymmetry production (for one flavor) per miit of phase space 
as a function of time. The red and blue lines correspond to the exact and perturbative results, 
respectively. Each plot covers a different time interval (time is expressed in GeV~^). We note that 
the exact and perturbative results agree very well up to t ~ 6000 ~ where we start to see 
small discrepancies that grow larger with time. 



frequency difference is produced by 1^12, ± — Ci;i2,±| 

m 

I {ml ± ml) — {mf ± mg) | 



\P\ 



MAM (i±^^ 



, (33) 

(e+5)K J 

which gives tg ~ 21000 GeV^^. Note that all estimates are roughly given by tg ~ \p\/MAM. 
As expected the two solutions follow each other nicely until t ~ 6000 ~ t^/S where 
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FIG. 4: Plot of the absolute error between the exact and perturbative results for two different time 
intervals. 

discrepancies appear. Those discrepancies are relatively small (roughly 10%), but they 
steadily grow with time; see Fig. (jl]) for plots of the absolute error between the two results. 
At t ~ ts ~ 21000 the absolute error becomes roughly as large as the absolute value of both 
asymmetries. We thus conclude that the exact and perturbative methods agree with each 
other from small times to a time tg at which secular terms are large and perturbation theory 
breaks down. 

V. CONCLUSION 

First principles computations of the lepton (baryon) asymmetry are difficult because some 
degrees of freedom need to be out-of-equilibrium in order to get a non-zero result, and the 
treatment of those degrees of freedom using quantum field theory is unwieldy. In this paper 
we have derived from first principles of quantum field theory and statistical mechanics a 
(simple) formula that could be used as a starting point for a perturbative computation of 
the baryon asymmetry of the universe. Our formalism is quite general and can be applied 
to other models. The only assumption entering into our derivation is that perturbation 
theory must be valid; physically this translates into the condition that the total number 
of out-of-equilibrium degrees of freedom (sterile neutrinos in our case) must remain small. 
This last condition depends on the parameters of the model under study. 

We have also tested this formula for the asymmetry on an exactly solvable toy model. 
We have confidence that the method works and that it can be applied to a more complicated 
model involving damping. The application of this formalism to the z/MSM and the study of 
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its phenomenology is work in progress. 
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